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Abstract 

In this work, we study a model of the chemostat where the species are present in two forms, isolated 
and aggregated individuals, such as attached bacteria or bacteria in flocks. We show that our general 
model contains a lot of models that were previously considered in the literature. Assuming that floccula- 
tion and deflocculation dynamics are fast compared to the growth of the species, we construct a reduced 
chemostat-like model in which both the growth functions and the apparent dilution rate depend on the 
density of the species. We also show that such a model involving monotonic growth rates may exhibit bi- 
stability, while it may occur in the classical chemostat model, but when the growth rate is non monotonic. 

Keywords. Chemostat, density dependent growth functions, flocculation. 

1 Introduction 

In culture of micro-organisms, the attachment of microbial individuals occurs frequently. The attachment 
can be either a "wall attachment" such as in the grovifth of biofilms or simply an aggregation such as in the 
formation of flocks or granules [S] |3H]- Flock or granule formation has a direct impact on growth dynamics, 
as the access to the substrate is limited for micro-organisms inside such structures. The mechanisms of 
attachment and detachment result from the coupling of hydrodynamics conditions and biological properties, 
but are not yet completely understood at the level of microbial individuals. Several attempts of computer 
models, using individual based representations, have been proposed and are under investigation for the 
simulation of these phenomenons, cf. for instance |16[ I34j . At a macroscopic level, substrate limitation 
can be measured experimentally in biofilms or flocks [H UHl SO]- A rough representation, suited to the 
macroscopic level, consists in splitting the overall biomass into two parts: a "planktonic biomass" , composed 
of free individuals and an "attached biomass" composed of individuals that are tied together [53]. This 
consideration leads to a signiflcant change on the performances predicted by the models, compared to purely 
planktonic cultures. In a chemostat-like device, planktonic cells are expected to consume easily the substrates 
necessary for their growth, but arc more keen to be carried out by the flow. On the contrary, cells among 
aggregates or biofilms have a more difficult access to the resources of the bulk fluid, but are more resistant 
to detachment induced by the hydrodynamical conditions. Therefore, mathematical models are expected to 
understand and predict the issues of these trade-offs. Several extensions of the well-known chemostat model 
[35j . considering two compartments of free and attached biomass for each species have been proposed and 
studied in the literature. In models with wall attachment, attached biomass is assumed to be fixed while 
detached individuals return directly to the planktonic compartment |321 136j . In models with aggregation. 
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aggregates are carried out by the flow but bacteria inside flocks are assumed to have no or reduced access 
to bulk resources [101 E] • 

Literature reports flocculation time scales of the order of 1 to 10 min [5121] to be compared with bacterial 
growth times of 1 h to 1 day, and with retention times of a few hours to a few days. Thus, considering that 
attachment and detachment processes may be fast compared to biological time, it is shown in [11] that the 
reduced dynamics of such systems amounts to have a single biomass compartment for each strain but with a 
density dependent growth rate. This justifies the consideration of density dependent growth functions in the 
chemostat model, as already introduced in the literature in the field of mathematical ecology [T] or waste- 
water process engineering [T2] . In [211 HSl [571 dH] , it has been shown that this could lead to the coexistence 
of several species in competition on a same limiting resource, thus invalidating the Competitive Exclusion 
Principle |13| (different mechanisms with considerations on the nutrient uptake could also lead to species 
coexistence [2]). 

In [11], the aggregates are assumed to have no biological growth (i.e. the attachment process is the only 
source of increase of the attached biomass). Aggregates are also assumed to be washed-out with the same 
dilution rate than planktonic cells. On the opposite, in wall attachment models, the attached biomass is not 
washed out at all. We believe that these two opposite cases (same dilution rate than planktonic biomass or 
no dilution rate) are too extreme to be fully realistic. In this paper, we revisit the chemostat model with 
two compartments, planktonic and aggregated biomass, but assuming that each biomass has its own growth 
rate and apparent dilution rate. This generalizes the two kind of models that we mentioned previously. 

The paper is organized as follows. In the Section 2, the general model of the system under interest is 
presented. In particular, for specific choices of different kinetics and mortality terms, it is shown that this 
model captures in fact many models of literature. In Section 3, this general model is reduced assuming the 
attachment/detachment processes are fast with respect to others. In the Section 4, the reduced model is 
analyzed for a specific class of models while its extension to the multi-species case and its analysis are carried 
out in Sections 5 and 6, respectively. Finally, discussion and conclusions are drawn in the last section. 



2 Modeling flocks or aggregates in the chemostat 

Consider the following model of the chemostat in which a population of microorganisms compete for a single 
growth- limiting substrate [HI |3T] : 

r ^ = D(5„ -S)- kf,{S)x 
\ X ^ {fi{S) - Dx)x 

In these equations, S{t) denotes the concentration of the substrate at time t; x{t) denotes the concentration 
of the population of microorganisms at time t] /i(-) represents the per-capita growth rate of the population 
and so y = 1/fc is the growth yield; Sin and D denote, respectively, the concentration of substrate in the 
feed bottle and the dilution rate of the chemostat; represents the removal rate of the population. 

Assume that the species is present in two forms: isolated or planktonic bacteria, of density u, and attached 
bacteria or flocks of bacteria, of density v. Isolated bacteria and flocks can stick together to form new flocks, 
with rate a{-)u, and flocks can split and liberate isolated bacteria, with rate /?(■)«: 

U > V, U i V. 

One obtains the following equations : 

r S = D{S,n~S)~f{S)u-g{S)v 

I u = {f{S)-Do)u-a{-)u + l3i-)v . (1) 
[ V = (giS) - Di)v + ai-)u - l3{-)v 

In these equations, S(t) denotes the concentration of the substrate at time t; u{t) and v{t) denote, respectively, 
the concentration of the population of planktonic microorganisms and flocks of bacteria at time t; f{S) and 



2 



g{S) represent, respectively, the per-capita growth rate of the populations; Sin and D denote, respectively, 
the concentration of substrate in the feed bottle and the dilution rate of the chemostat; Dq and Di represent, 
respectively, the removal rate of the microorganisms. 

The dot in attachment rate a(-) and detachment rate means that these rates can depend on the state 
variables, so that system ([T|) recovers some of the models which were considered in the existing literature. 
For instance, the model of adaptive nutrient uptake, where u denotes the low growing cells and v denotes 
the fast growing cells considered in [37j is obtained with attachment and detachment rates depending only 
on S 

a(.)=a(5), /?(•)= 

The model of wall growth, where u denotes the density of planktonic bacteria, and v denotes the density of 
wall-attached bacteria, considered by Pilyugin and Waltman [32|, is obtained with constant rates 

a{-)=a, Pi-)=b. 

The Freter model |S] is given by 

S = D{S„,-S)- f{S)u-g{S)v 

u = {f{S)-Do)u-a{l-W)u + bv + g{S){l~G{W))v (2) 
V = {g{S)G{W)-Di-b)v + a{l-W)u 

where W = w/wmax and G{-) is decreasing. Notice that this model is a particular case of the model ([T]) with 

«(•) = a(l - W), /3(.) = b + - G{W)) . 

Actually, if Wmax = oo one obtains W = and if G{0) = 1 then a(-) ~ a, f3{-) ~ b, and is simply the 
model of Pilyugin and Waltman 02] • 

A model with flocks of two bacteria has been considered in [11] 

S = D{S^n-S)- f{S)u~g{S)v 

ii = {f{S) - Do)u-au^ + bv (3) 
i, = [g{S) ~ Di)v + av? - bv 

that is a particular case of model ([ij obtained with 

a(-) = au, I3{-) ~ b. 

This model has been studied by Haegeman and Rapaport [TI] in the case of g{-) = where the bacteria in 
flocks are assumed to do not consume any substrate, and by Fekih-Salem and al [7] in the more general case 
ofO^(7(-)^/(-) where the bacteria in flocks consume less substrate than the isolated bacteria. This model 
has been also extended to the case of flocks with an arbitrary numbers of bacteria in [10] . 

In the present paper we will not consider the size or the number of bacteria in flocks in our model. We 
simply distinguish the biomass in flocks and the isolated biomass. The biomass of isolated bacteria is denoted 
by u and the biomass in flocks is denoted by v. Hence isolated bacteria and isolated bacteria or flocks can 
stick together to form new flocks, with rate a(u-\-v)u, where a is a constant, proportional to both the density 
of isolated bacteria, that is u, and the total biomass density, that is u + w, and flocks can split and liberate 
isolated bacteria, with rate bv, where 6 is a constant, proportional the their density v. Hence, taking 

«(•) = a(u + u), P{-)^b 

in model ([T]), one obtains the following dynamical system 

S = D{S^n-S)- J{S)u-g{S)v 

u = ifiS) - Da)u-a{u + v)u + bv (4) 
ii = (giS) — Di)v + a{u + v)u — bv 

that we aim to study in the present paper, as well as its extensions to multi-species populations. 
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3 A two time scales dynamics 

The general model for the floeculation is 

D{S,^ -S)^ f{S)u - g{S)v 
{f{S) — DQ)u — a{S,v,u)u + l3{S,u,v)v . 
{g{S) — Di)v + a{S, v, u)u — l3{S, u, v)v 

If we assume that the dynamics of floeculation and defloeculation is much faster than the growth of the 
species, one can write the model the following way 

f S = D{S,n-S)-f{S)u-g{S)v 

ffrQ\ n \ a{S,u,v) /3{S,u,v) 
u = {,f[S)-Dn)u itH V (5) 

V = (g[S) — Ui)v-\ u V 

e e 

where e is expected to be a small no n- negative number. Notice that the dynamics of the total biomass 
a; = u + u is given by the equation 

x = {f{S)-D^)u+{g{S)-D^)v. 
Thus, u and v are fast variables, while S and x are slow ones. The fast dynamics is given by 

u' ~ —a(S, u, v)u + I3{S, u, v)v , 
and the slow manifold is defined by the positive solutions of the system 

a{S^ u, v)u — I3{S, u, v)v with u + v = x . 

Hence one has 

u=p{S,x)x, V = {1 — p{S,x))x. 

Assuming that this slow manifold is asymptotically stable for the fast equation, the reduction of the system 
to the slow system gives an approximation of the solutions of (O for arbitrarily small e, in accordance with 
the theories of singular perturbations or variables aggregation 30, 3, 33^. The reduced model is obtained by 
replacing the fast variables u and v in the equations of S and x: 

S = D{S,n - S) - ^(5, x)x 
X = {ij{S, x) — d{S, x)) X 

where 

f,{S, x) = x)f{S) + (1 - p{S, x))giS), 

d{S, x) = piS, x)Do + (1 - p{S, x))Di. 

Notice that p(-) depends on functions a{.) and /?(.). Consequently the growth function fi{-) is a density 
dependent growth function model, as already studied in |26 1 [27 1 128 ] . But, here the removal rates are replaced 
by a function d{-) that depends also on functions a{.) and /?(.). This last property is new in the literature, 
to our knowledge. 

For the slow and fast case of ([3]) one has, see [11] 

^ = ^ and pix)^ ' 



e £ e e 1 + + 4fx 



4 



The slow and fast case of ^ is given by 



S = DiS,n-S)-f{S)u~g{S)v 

u={f{S)-Do)u-^(u + v)u+!^v . (7) 
V = {g{S) - Di)v + f (u + v)u - 



The slow manifold is defined by the positive solutions of the system 

a{u + v)u = bv with u + v = x . 

Hence one has 

u = p(x)x, V = (1 — p(x))x, where v(x) — . 

Hence, the reduced model corresponding to is given by 

S = D{Sin ~ S)- fi{S, x)x 
X = x) — d(x)) X 



(8) 



where 

lo \ bf{S) + axg{S) hD^ + axDi 

fi{S,x)^ , d{x)^ — . (9 

+ ax + ax 

In the rest of the section, we use Tikhonov's theory [29l |39l |42] (see also [20]), and we justify that the 
solutions of ([7]) are approximated by the solutions of the reduced model (|8l9p . 

Theorem 3.1 Let (^(i, e), e), e)) be the solution of ^ with initial condition (S'o,uoi^'o) satisfying 
So ^ 0, uq > 0, and vq ^ 0. Let {S{t),x(t)) be the solution of the reduced problem with initial conditions 

S{0) = So, x{t) ^uo + vq. 

Then as e — > 0, 

S{t,e) = Sit)+o{l), x{t,e)^x{t)+o{l) (10) 

uniformly for t G [0, T], and 

b ax(t) 
u{t,e)^——-x{t) + o{l), v{t,e)^ x(t) + o(l) (11) 

b + ax(t) b + axit) 

uniformly for t G [to,T], where < to < T are arbitrary but fixed and independent of e. If the solution 
{S{t),x(t)) of the reduced problem converges to an asymptotically stable equilibrium, then we can put T = +oo 
in the the approximations ilO\} and ill]} given. 

Proof. In the variables S, x ~ u + v and u, system ([7]) is written 

S = D{S,n ~S)- f{S)u ~ g{S){x - u) 

X = f{S)u + g(S)(x — u) — Dqu — Di(x — u) /-i 2~) 

ah 

u = {f{S) — Do)u xu -\ — {x — u) 

This is a slow and fast system with two slow variables S and x and one fast variable, u. The fast equation is 

u' = —axu + b{^x — u) = hx — {ax + b)u. (13) 
The slow manifold (or quasi-steady-state) is given by 

b 



u = xp{x), where p{x) 



ax 
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Since ax + b > this slow manifold is globally asymptotically stable. Thus, the Tikhonov's Theorem 
[29l [39l |42] applies and asserts that after a fast transition toward the slow manifold, the solutions are 
approximated by a solution of the reduced equation, which is obtained by replacing the fast variable u in 
([T^ by the quasi steady state u = xp{x). This reduced system is 

S = D{S,n -S)- {f{S)pix) + g(S){l - p{x))) x 
X ^ if{S)pix) + .9(5)(1 - pix)) - Dop{x) -D,{1~ p{x))) X ' 

This is system (|8l9p . The approximations (|10p and pT|) follow from the Tikhonov's Theorem. Recall that 
when the reduced problem has an asymptotically stable equilibrium, then these approximations hold for all 
t > and not only on a compact interval [0,T). Recall also that there is a boundary layer for the fast 
variables u and v, that is the approximations hold only for t ^ to where to > can be arbitrarily small 
but fixed. Therefore, very quickly, the density v{t) of flocks tends to 

a(x(0))2 



b + ax{0) ' 

and the density u{t) of planktonic bacteria tends toward 

6a;(0) 
b + ax{Q) ■ 

These values depend only on the initial total biomass x{0) — u{0) + f (0), and on the ratio a/b between the 
rate of flocculation and deflocculation. After this boundary layer, one has a slow variation of the densities 
of flocks and isolated bacteria according to the density dependent reduced model (|8l9p . 

4 Study of the reduced model for one species 

We consider in this section the mathematical analysis of the reduced model ([8]) . We do not assume that the 
functions fi{S,x) and d{x) are of the particular form We assume that 

HO: /i(0, x) = and yu(5, x) > for aU S* > and all x ^ 0. 

HI: 1^ > and ^ < for all > and all x ^ 0. 

dS ox 

H2: rf(0) = Do, d{+oo) = Di < L», d(x) > 0, d'{x) < and [xd{x)]' > for ah x ^ 0. 

Let us denote by 

fiS) = 0), and giS) - ^l{S, +oo). 

The functions /(•) and g{-) arc increasing and positive for all S > 0. If equations f{S) = Do and g{S) = Di 
have solutions, one lets 

Ao = .r'(^o), and Xi = g-\Di). 
Otherwise one lets Xk = +oo, fc = 0, 1. We add the following assumption 

du 

H3: If Ao < Ai, then for all S £ [Aq, Ai) and x ^ one has d'(x) > -^(S*, x). 

ox 

dii 

H4: If Ai < Ao, then for all S e (Ai,Ao] and .t ^ one has d'(x) < ■^(S,x). 

ox 

Lemma 4.1 Assumptions H0-H4 are satisfied in the case 

^i{S, x) = f{S)p{x) + g{S){l - p{x)), d{x) = Dop{x) + D,{1 - p{x)) 

where /(•) > g{-) are increasing functions and p{x) is a decreasing functions such that p{0) = 1, p{+oo) ~ 
with [xp{x)y > 0. 
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Proof. Since p'{x) < and f{S) > g{S) for all S > 0, one has 



dx 



{f{S)~g{S))p'{x)<0. 



Thus HI is satisfied. On the other hand 

d'{x) = {Do - Di)p'{x) < 0, [xd{x)]' = Di + (L>o - Di)[xp{x)]' > 
since p'{x) < 0, [xp{x)]' > and Dq > Di. Thus H2 is satisfied. Moreover 

- d'ix) = [f{S) -Do + D,- g{S)]p'{x) < 

for ah S e [Ao, Ai) and x ^ 0, since p'{x) < and f{S) ^ Dq, Di > g{S) for Aq < 5 < Ai, sec Figure [1] 
Thus Assumption H3 is also satisfied. Similarly one has 





Ao Ai Ai A() 

Figure 1: On the left, the case Aq < Ai. On the right, the case Ai < Ag. 



dx 



{S,x) - d'{x) = [f{S) -Do + D,- g(S)]p'{x) > 



for ah S e (Ai, Ao] and x ^ 0, since p'{x) < and f{S) ^ Dq, Di < g{S) for Ai < 5 ^ Ao, see Figure [H 
Thus Assumption H4 is satisfied. 



Notice that, ii p{x) 



b+ax 



then the properties of p(-) stated in the lemma are satisfied. Indeed 



p'{x) = 



—ah 



[h + axY' 



[xp{x)]' = 



[h + axY ■ 



4.1 Existence of equilibria 

The equilibria of the system are solutions of the set of equations 

D{Sm - S) - fj.{S, x)x = 
{fJ,{S, x) — d(x))x — 

The second equation is equivalent to x = or x) — d{x). If x = then from the first equation one has 
'S' = Sin- This is the washout equilibrium 

Eq = (S',;„, 0). 
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If i-i{S, x) = d{x), the first equation gives D(Sin ^ S) ~ xd(x). Hence 

xd(x) 



S = j{x) := S^., 



D 



Since 



7(0) = S,„, and j'{x) = < o 



the function 7(-) is decreasing. Thus one have to solve the equation 

/i(S', x) — d{x). 

Since ^ > 0, by the implicit function theorem, thus equation defines a function 

5* = such that Ao — 0(0), 

and 



(14) 



c^'ix) 



d'{x) - |^((/)(x),a;) 



The sign of </>'(•) is given by assumptions HI, H3 and H4. The cases Ao < Ai and Ao > Ai have to be 
distinguished. 

When Ao < Ai the function S ~ (/)(x) is defined for all x ^ and satisfies 

Ao = </>(0), Ai==0(+oo), </.'(x)>0. 
The equilibria are the intersection points of the graphs of functions 

S = (j){x) and S ~ ^[x). 

Since the first function is increasing and the second one is decreasing, there is a unique solution if Ao < Sim 
and no solution if Ao > Sin , see Figure [2] 



S = (t){x) 




S = j{x) 



Ao Sin 



Sin Ao 



S 



Figure 2: NuU-clines S = (f){x) and S = 7(0;) in the case Aq < Ai. On the left, the case Ao < Sin with a 
unique intersection point. On the right, the case Ao > 5*^,1 with no intersection point. 

Proposition 4.2 //Ao < min(Ai, S'm), there exists a unique positive equilibrium. If Sin < Ao < Ai, there is 
no positive equilibrium. 

When Ai < Ao the function S = Mx) is defined for all x ^ and satisfies 



Ao = <^(0), Ai=(/.(+w), c^'{x)<Q. 
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Figure 3: Null-clines S = (j){x) and S = 7(2;) in the case Xq > Xi. On the left, the case Xq < Sin with at 
least one intersection point. On the right, the case Aq > Sin with no intersection point. In the second case 
it is possible to have two intersection points, see Figure 21 

Both functions S = 4>{x) and 5* — ^{x) are decreasing and 

0(0) Ao, S,n = 7(0). 

Thus, if Aq < Sin the graphs intersect at at least one non- negative point (sec Figure |3]). If Aq > Sim the 
graphs of functions S = (t>{x) and S = 7(x) can interest or not (see Figures |3] and |4]) . 

Proposition 4.3 //Ai < Aq < Sin, then there exists at least one positive equilibrium. Generically one has 
an odd number of positive equilibria. If Xi < Sin < Aq , then the system has generically no positive equilibrium 
or an even number of positive equilibria. 




Figure 4: On the left, the growth functions /(•) and g{-) given by (|15|) . On the right, the corresponding null- 
clines S = 4>{x) and S — ^{x) for the parameter valucs (jl6p . showing the existence of two positive equilibria 
in the case when Aq > Sin. 



4.2 Stability of equilibria 

The .lacobian matrix of ([5]) is given by 

—u — X 

J 



x^{S,x) 



~a;|^(S',a;) - ^(S", x) 



x^{S,x) 



^j.{S,x) 



dx 



d{x) + .T— ^(5, x) — xd'{x) 
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At washout Eq = (S'i„,0) this matrix is 

r -D -f{S,n) 

- Do J 

The eigenvalues are —D and f{Sin) — Dq. Hence one has the following result 

Proposition 4.4 If Xq < Sin, then Eq is unstable (saddle point). If Xq > Sin, then Eq is locally exponen- 
tially stable (stable node). 

At a positive equilibrium Ei = {S,x), one has necessarily fj.{S,x) = d{x), and the Jacobian matrix is 



Jl 



-D-Xj^{S,x) -x^{S,x) - fi{S,x) 



x/g{S,x) 



x-^{S,x) —xd'{x) 



The trace of Ji is 



Since = d{x) + xd'{x), one has 

tr(Ji) = -D + d{x) ™ x^{S,x) + x?l^{S, x) - < 

ab ox 

since d{x) < £), || > 0, §f < and [xd{x)]' > 0. The determinant of Ji is 



det(Ji) = Dx (^d'{x) + 



One can write this determinant as 



det(Ji) = Dx-^ 



5m , [xdjx)]' 



dS 



D 



^Dx^i<P'{x)-j'ix)). 



Thus, If (f>'{x) > 7'(a;), then the determinant is positive and hence the eigenvalues are of negative real parts: 
the equilibrium Ei is locally asymptotically stable. On the other hand, if (j)'{x) < ^'{x), then the determinant 
is negative, and the eigenvalues are of opposite sign: the positive equilibrium is a saddle point. Thus we 
have shown the following result 



Proposition 4.5 The following cases occur. 

1. If Xq < min(Ai, S',;„), the unique positive equilibrium is locally asymptotically stable (and fulfills 4>'{x) > 
> i{x)). 

2. If Xi < Xq there are possibilities of multiple positive equilibria, that are saddle points when 4>'(x) < ^'{x) 
or stable nodes when 4>'{x) > j'{x): 

- when Xq < Sin, there exists at least one positive equilibrium and one has an odd number of 
equilibria which are alternatively stable and unstable, 

- when Xq > Sin, there is no or an even number of equilibria which are alternatively stable and 
unstable. 
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One concludes that in the case Ai < Aq, the system can exhibits bi-stabihty of the washout equilibrium 
and a positive equihbrium. If the initial density of flocks is small enough, the solutions will converge to the 
washout equilibrium, otherwise it will converge to the positive equilibrium, sec Figure [5l In this figure we 
notice that the washout equilibrium is stable together with the positive equilibrium corresponding the lowest 
value of 5*. The domains of attraction of the stable equilibria arc separated by the stable separatrix of the 
positive saddle node. The simulations shown in Figures S] and [5] where obtained for (|8l9p . with the following 
Monod functions 

,(S)^J^ (15) 

and the following values of the parameters 

Do = D = l, Di = 0.5, a = 4, b=l, S,,, = 0.9. (16) 




Figure 5: For the Monod functions ([T5|) and the parameters values the system exhibits bi-stability. 



5 Flocculation with several species 

We assume that n species are competing on a same limiting resource, and that each species is present in 
two forms: isolated bacteria, of density Ui, and bacteria in flocks, of density Vi, for i — 1 ■ ■ ■ n. We assume 
that isolated bacteria can stick with isolated bacteria with flocks to form new flocks, with rate ai{-)ui. We 
assume also that flocks can split and liberate isolated bacteria with rate j3i{-)vi. 

Ui > Vi, Ui < Vi . 

Then the equations are 

n 

S = Z?(S„ - 5) - Y^iMSy^ + g^iS)v,) 

1=1 (17) 
Ui = {MS) - Doi)ui - ai{-)ui + l3i{-)vi, 1 ^ i ^ n 

V, = {gi{S) - Dii)vi + ai{-)ui - (ii{-)vi 
The dynamics of the total biomass densities Xi ~ Ui + Vi of the species i is 

Xi = fi{S)ui + gi{S)vi - DoiUi - DuVi . 

We consider here the case where 

n 

ai{-) = ^AijXj, Pi{-) = Bi 
where Aij and Bi are non- negative constants. 
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Remark 5.1 By letting the Junctions oii{-) depending on Xj with j ^ i, we implicitly consider that flocks or 
aggregates can incorporate individuals of different species. 

The removal rate Du has to be interpreted as the sum of the removal rate of the aggregates and the mortality 
rate of each species i in its attached form. 

Moreover, we assume that the dynamics of flocculation and dcflocculation are fast compared with the 
dynamics of the growth of bacteria, that is 

A - ^ B - — 
In the variables 5, x^, Ui, the system is written 

n 

S = D{S,n - S) - + gi{S){xi - Ui)] 

1=1 

ii = fi{S)u.j,+ g^{S){x^-u.i)- DoiUi- Du{x^-Ui), i 1 ■ • • n 

1 " b 
= {fi{S) - Doi)ui ~ 7 ai]XjUi + —{xj - Ui) 



(18) 



This is a slow/fast system with the variables S*, Xi are slow and the variables itj fast. The fast equations are 

n 

u[ = — ''^^OijXjUi -\- biixi — Ui)^ i = l---n (19) 
where Xi are considered as parameters. The slow manifold (or quasi steady-state) is given by 

biXi 



1 ■ ■ ■ n . 



(20) 



bi H~ ^ ^ dijXj 

Since one has X]J=i '^ij^j 0, for i = 1 ■ ■ ■ n, this slow manifold is globally asymptotically stable for ^ 
Thus, the Tikhonov's Theorem [29l [39l |42] applies and asserts that, after a fast transition toward the slow 
manifold, the solutions are approximated by a solution of the reduced equation, which is obtained by replacing 
the fast variables Ui in (fT8|) by the quasi steady states ([20|) . One obtains the following reduced model 



(21) 

(22) 
(23) 



S = D{Sin ^ S) ^J,i{S, x)xt 

1=1 

x.^ = {fii{S,x) - di{x)) x^, i = l---n 

where 

fii{S,x) = fi{S)pi{x) + gi{S){l -pi{x)), di{x) = DoiPi{x) + Du{l ~pi{x)), i = 1, • ■ • ,n 
with ^ 

Pii^) = ; X = {Xi,- ■ ■ ,Xn)- 

bi + OijXj 
Hence, we have shown the following result 

Theorem 5.2 Let {S{t,e),ui{t,e),vi{t,e), ■ ■ ■ ,Un(t,e),Vn{t,e)) be the solution of |_?7| ) with initial condi- 
tions S{0) ^ and Ui{0) > 0, Vi{0) ^ 0, for 1 ^ i ^ n. Let (^S{t),xi{t), ■ ■ ■ ,x„(t)) be the solution of the 
reduced problem i21\) with initial conditions 



S{{)) = 5(0), x^{t) = u,{Q) + v,{{)), l^i^ 



12 



Then as £ — !• 

S{t,e)='S{t)+o{l), Xi{t,e)=x,{t)+o{l), l^i^n 
uniformly for t G [0, T], and for all 1 ^ i ^ n, as £ — > 

biX^(t) (Z]"=i 

Ui{t,£) = , I v^n ' ^777+0(1), v^{t,e) = ^ „ zr77^ + o(l) 

uniformly for t € [to,?^]; where T > > are arbitrarily fixed. If the solution of the reduced problem tends 
to an asymptotically stable equilibrium, then we can put T = +00 in the approximations given above. 

Since the planktonic bacteria have a better access to the substrate than the bacteria in flocks one assumes 
fi{S) > gi{S). Notice that one has ^ < for any i, j. Hence 

|l.,/.(S,-,,5,)|i<0 

with 

^ = fl{S)p,{x) + gl{S){l - p,{x)) > . 

As for the one species case this approach give a motivation to density dependent growth function models, 
that may lead to species coexistence [26l [27l [28] . 

6 Study of the reduced model for several species 

The results of this section apply to the reduced model (|21I22I23|) in the particular case where the rate of 
attachment and detachment of species Xi with species Xj are negligible for i ^ that is to say Uij = for 
i ^ j ■ In that case the function pi depends only on Xi and is given by 

/ \ ^» 
Pi[x) 7— 

SO that the growth function fii and removal rates di in (|2ip depend only on Xi. We consider then the model 

n 

S = D{S.in - S) ~^fi.i{S,x^)x^ ^24) 

Xi ^ [^J.i{S, Xi) - di{xi)]xi i = 1, • • • , n 

This model was studied in |26j . in the case when di{xi) — D. We do not assume that the functions fj,i{S,Xi) 
and di{xi) are of the particular form given by (|22p . We assume that 

H5: /ii(0, Xi) = and /ii(S', Xi) ^ for all S* > and all Xi ^ 0. 

H6: ^ > and ^ < for aU S* > and all x, ^ 0. 

Ob oxi 

H7: di{Q) = Dai, di{+oo) = Du < Dm ^ D, di{xi) > 0, d[{xi) < and [xidi{xi)]' > for aU Xi > 0. 
Let us denote by 

fi{S) ^ fii{S,0) and gi{S) ^ fii{S, +00). 

The functions fi{.) and gi{.) are increasing and positive for all S > 0. If equations fi{S) — D^i and 
gi{S) ~ Dii have solutions, one let 

Aoi = f^i^^iDot) and Ah = g^^iD^) 

otherwise one let \ki = +c«, fc = 0,1. As for the case of one species (see Assumption H3), we add the 
following assumption 
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H8: Aoi < Xii for i = 1 • • - n. For all S G]Aoi, Aii[ and Xi ^ 0, one has d^{xi) > ^{S,Xi). 

If the inequality Aoi < Aii is reversed for some i ~ 1 ■ ■ ■ n, then the situation is mueh more diffieult and will 
be studied in the future. Denote 

Ao = max{Aoi; i = 1, - ■ ■ , n} and Ai = min{Aii; i = 1, - ■ ■ , n}. 

We assume that 

H9: Ao < min(Ai, S',;„). 



We consider here the existence of a positive equilibrium. The equilibria of (|24p are solutions of the set of 
equations 

/ n 

D{Sin ~ »S') ~ ^ ^ fJ'i{S, Xi)Xi 

i=l ^^^> 

^ii{S,Xi) = di{xi) or Xi=Q i=l,---,n. 
Thus we have to solve the equations 

^j.i{S,Xi) = di{xi). 

Since H6, by the implicit function theorem, this equation gives a function S = 4>i{xi) defined for all Xi ^ 0, 
such that 0i(O) = Aoi, (/)i(+oo) = Xu and 

= Qj:^ > 0. 

The sign of (/>'{■) is given by assumptions H8. We define the function Xi : S i-^ Xi{S) on [0, Xii[ by 
X^:[0, Xi,[ M+ 

' =0,ri(5) if Ao, ^^<An. 

Let hi{S) — ^ii{S, Xi{S))Xi{S). Since Xi{-) is increasing over [Aoi, Aii[, so is hi{-) over this interval. Indeed, 
one has 

KiS) ^(^ + ?!^X^{S)) X,iS) + {S,X,{S))XliS). 



^ dS dx 

Moreover, for S e]Ao„Ai,[, ^i, iS,X,{S)) - and 

X^{S) = ^ ^ > 0. 

4(x,{s))~£-is,x,is)) 

Then 

h',{S) = {X,{S)) X,{S) + {X,{S))] XiiS), for S e]Ao„ Ai,[. 

Using H7, 

[d.i{xi)xi]' d[{xi)xi + di{xi) > 0, for Xi ^ 0. 
Hence the sign of ft.-(S') is the same as the sign of -'^^,'(5'), that is, hi{-) is increasing over [Ao;, Aii[ (see Fig. 

Consider now the function 



=^/i,(5)-i?(5,;„-5). 



1=1 
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Figure 6: Condition of existence of the positive equilibrium of p4| for n = 3. 
Lemma 6.1 Equation H{S) = admits a unique solution S* G]0, Ai[. 

Proof. Since hi{S) = for 5* G [0, Aoi] and hi{S) is increasing over [Aoi, Xii[ (see Fig. the function H{-) 
is increasing over (0,Ai), and 

H{0) = -DS^n < and lim H{S) = +oo. 

Hence, there exists a unique S* g]0, Ai[ such that H{S*) = 0. We have the fohowing resuh : 
Proposition 6.2 Assume that H5-H9 hold. System \24-^ has a unique positive equilibrium if and only if 

n 

^Ai,(Ao,X,(Ao))x,(Ao) < D{S,r. - Xo)- (26) 

i=l 

Proof. A positive equihbrium E* = {S*,xl, ■ ■ ■ must satisfy 

n 

I?(5™-^*) = ^//,(^*,a;*K (27) 

and 

/i.(5*,<) = d.(x*). (28) 
Equation (f^ is equivalent to = Xi{S^). Thus, (P7|l can be written 

n n 

1=1 1=1 
that is ^^(5"*) = 0. Since J2i=i htiS*) > 0, then one must have 

S* < Sin and S* > Xq. 

Notice that Aq < S** < Ai and S* < Sin arc satisfied if H9 holds. Then, since ^(5") est increasing over 
[0,Ai[, 

Xo<S* ^ hCXo) < H{S*) = 0. 
Therefore there exists a unique positive equilibrium S* exactly when H{Xo) < 0, which is equivalent to (j26p . 

We study now the asymptotic behavior of the positive equilibrium. First, we establish the following result 
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Lemma 6.3 Consider the matrix 





-D - J27=i 


Cl 


C2 






ai 


-bi 


• 





A = 


ai 





-b2 ■ 







an 





• 


■ • -bn _ 


and for i = 1 ■ 


■ n, Oi ^ 0, bi > 


and 


Cl s^bi. 


Then all 



(29) 



real part. 

Proof. Let A be an eigenvalue of A and V = {vq^vi, - ■ ■ ,Vn) ^ a. corresponding eigenvector. Hence we 
have 



-D - ^ a, vq + ^ c,u, = Xvq 



(30) 



aiVQ ~ biVi — Xvi i — 1, ■ ■ ■ ,n. 
Assume that a = Re(A) ^ 0. Since bi > 0, then X + bi 0. Therefore, 

for i = I ■ ■ ■ n. 



A + b, 



If Wo = then u,; = for i = 1 • • • n, so that V ^ which is impossible. Thus vo ^ and from the first 
equation of pOI) . we deduce, after simplification by vq that 



i=l 4=1 

Let X = a + ifi. Taking real part of both sides one obtains 



A + 6, 



A. 



Ci{bi + a) 
f'{b~Taj^TW^ 



Since Ci ^ bi and 6; + a > then 



Ci{bi + a) ^ bi{b, + a) ^ bjjh, + a) _ b, ^ ^ 



Since a,- ^ then 



Hence 



(&, + a)2 + /?2 [bi + a)2 + /?2 + a)2 6, + a 

Ci{bi + a) r ■ -, 

n n 

a ^ -D - E + < 

i=l i=l 



which contradicts a ^ 0. 



Then, we state the following result : 
Proposition 6.4 // E* exists, then it is locally exponentially stable. 

Proof. Since fii{S*,x*) — di{x*), the Jacobian of the system ([M)) at E* is of the form where 

= "^^('5'*, 2;*)a;* 6; = — ^— (5*, a;*).Tj + a:*c?^(x*), Ci = —-^—{S' ,Xi)xi ~ di{xi). 

Since H6, a.^ > 0. Since H7, (ii(x*) + a;*d-(x*) > 0, then -di{x*) < x*d^{x*) and hence Ci < bi. Since H8, 
bi > 0. The result follows from Lemma [6.31 
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7 Discussion and conclusion 



In this work, wc considered a general model of a bio-process with three compartments involving the sub- 
strate, the planktonic and the attached biomass densities, respectively. Each compartment of the biomass 
is characterized by its own specific growth rate and apparent dilution rate, generalizing previous models of 
biofilms (with no dilution rate for the attached bacteria) or models of perfect flocks (with no growth rate for 
aggregated individuals). We have analyzed a class of such models with planktonic and structured biomass, 
under the assumption that attachment and detachment processes are fast compared to the biological scale. 
Notice that it is only under this assumption that the main results of the paper are valid, notably the fact 
that the three order model ([T]) can be reduced to the second order model ©■ If it is not the case, the analysis 
of the original three order model must be done to establish its qualitative behavior. 
Our study reveals two main characteristics of this model: 

1. the reduced dynamics may exhibit a bi-stable behavior even though each growth function is monotonic. 
This phenomenon is new and is usually met in the chemostat but when the growth function is non- 
monotonic (such as the Haldane law); 

2. for bio-processes in which part of the biomass is under a structured form (in flocks or biofllm), the 
macroscopic models (with reduced dynamics involving only the aggregated biomass and the substrate) 
should include a growth rate and an apparent dilution rate that are both density dependent. This 
result contributes to the actual debate in biotechnological engineering involving bio-processes with 
structured biomass, where it was not clear whether it is better to modify the growth rate functions or 
the hydrodynamical terms in the macroscopic equations of the system, to cope with the speciflcity due 
to the attachment process. 
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